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Abstract 

It is shown that the a 2 — dynamo of Magnetohydrodynamics, the hydrodynamic Squire equation as well 
as an interpolation model of VT— symmetric Quantum Mechanics are closely related as spectral problems in 
Krein spaces. For the a 2 — dynamo and the VT— symmetric model the strong similarities are demonstrated 
with the help of a 2 x 2 operator matrix representation, whereas the Squire equation is re-interpreted as 
a rescaled and Wick-rotated VT — symmetric problem. Based on recent results on the Squire equation the 
spectrum of the VT— symmetric interpolation model is analyzed in detail and the Herbst limit is described 
as spectral singularity. 

1 Introduction 

Non-Hermitian VT— symmetric quantum mechanical systems 0I20SE1I11IZ1IHI are known to possess spectral 
sectors with purely real eigenvalues as well as sectors with pairs of complex conjugate eigenvalues. Changes of 
certain system parameters can lead to spectral phase transitions from one sector to the other. The physics in the 
two sectors has been identified with phases of unbroken VT— symmetry (real eigenvalues) and spontaneously 
broken VT— symmetry (pairwise complex conjugate eigenvalues) ^|2]. From a mathematical point of view, 
non-Hermitian VT— symmetric Hamiltonians are self-adjoint operators in Krein spaces IIUI III] — Hilbert 
spaces with an additional indefinite metric structure — and the two spectral sectors correspond to Krein space 
states of positive or negative type (real eigenvalues) and neutral (isotropic) states (pairwise complex conjugate 
eigenvalues) . 

Apart from VT— symmetric Quantum Mechanics (PTSQM), it is known that a certain class of spherically 
symmetric mean-field dynamo models |12| of Magnetohydrodynamics (MHD) can be described by self-adjoint 
operators in Krein spaces as wel l |13| . These models show similar spectral phase transitions from real to pairwise 
complex conjugate eigenvalues |14| — and only the physical interpretation differs from that in PTSQM. For 
dynamos it simply consists in a transition from non-oscillatory states to oscillatory states. 

In the present paper, we are going to briefly describe the underlying structural operator theoretic parallels 
between PTSQM models and the spherically symmetric MHD a 2 — dynamo (section The discussion will be 
illustrated with the help of a VT— symmetric interpolation between a harmonic oscillator placed in a square well 
and an empty square well (section [3J . This interpolation shows a rich structure of spectral phase transitions 
with a couple of unexpected features. Furthermore, we will show in sect. 0]that the eigenvalue problem of 
the VT— symmetric (intermediate) interpolation model with linear complex potential (purely complex electrical 
field) within the square well is mathematically identical to the eigenvalue problem of the rescaled and Wick- 
rotated Squire equation of hydrodynamics which describes the normal vorticity of a plane channel flow (Couette 
flow) with linear transversal velocity profile. Recent Airy function based results on the Squire equation allow 
us to analytically describe the spectral behavior of the PTSQM model in the limiting case when the width of 
the square well tends to infinity. In this limit we reproduce the Herbst model Xh with its empty spectrum. 
The limiting behavior occurs as a blowing-up of the spectrum to infinity along three directions on the complex 
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plane — leaving behind a spectrally empty region at any fixed finite distance from the origin of the spectral 
plane. In section [3] we briefly sketch some links of the obtained results to other physical setups and analytical 
techniques. 



2 Krein space properties of VT— symmetric quantum models and of 
the spherically symmetric MHD a 2 — dynamo 

2.1 VT — symmetric quantum models 

In their seminal letter Bender and Boettcher identified VT— symmetry as the essential property of the 
non-Hcrmitian quantum system 

d 2 

Hip(x) = Eip(x), H = — — +gx 2 (ixY (1) 

ax* 

which ensures the reality of its spectrum for exponents v <E [0, 2) and ip(x) € H = L2(—oo, oo) [TS]. This allowed 
them not only to extend an earlier conjecture of Bessis and Zinn- Justin (whose numerical results indicated that 
quantum systems with complex potential V(x) = ix 3 might have a purely real spectrum), but also initiated 
the still lasting intensive study of generalized VT— symmetric non-Hermitian systems |17| . Such systems are 
characterized by a VT— symmetric Hamiltonian H, 

[VT, H]=0 (2) 

where V denotes a reflection 

VxV = -x, Vip(x) = ip(-x) (3) 
while the time-reversal operator T performs complex conjugation 

TiT = -i, Txjj{x) = xjj{x)*. (4) 
Because both operators V and T are involution operators, 

V 2 =1, T 2 = I, (5) 

they induce natural Z2— gradings of the Hilbert space Ti. For our subsequent analysis it suffices to consider 
the subclass of models which can be defined solely over the real line a; 6 ML For such models the T— induced 
Z2— grading corresponds to a splitting of the wave functions ip G TL into real and imaginary components (what 
is of no direct physical interest in a quantum mechanical context; additionally one would have to work in a real 
Hilbert space with doubled dimension compared to the original complex one), whereas V induces a Z2— grading 
into parity even and parity odd components 

ip(x) = ip+(x) + ip-(x), Vip±{x) — tp±(— x) = ±ip±(x) . (6) 

The corresponding Z2— graded Hilbert space splits as 

H = H+®H-, V±eW±. (7) 
In the case of a simple VT— symmetric one-particle system with Hamiltonian 

H = -dl + V+{x) + iV-(x), V±(-x) = ±V±(x), ZV± = (8) 

it holds 

H = VH f V (9) 

and V is a so called fundamental (canonical) operator symmetry 0^]] of H — i.e., H is V— pseudo-Hermitian 
in the sense of Refs. Operators with an involutive fundamental symmetry are known to be symmetric 

— and for appropriately chosen domain (boundary conditions for the functions ip{x)) even self-adjoint — in a 
Krein space JC. For V— pseudo-Hermitian operators over the real line this Krein space IC-p is given as jl II 1181 [TT?| 



(/C P , [.,.]„), [^^} v = (ip,V<t>)= i>*(x)Vcl>(x)dx= r{x)^{-x)dx. (10) 

JCCR JCCR 

Depending on the concrete problem, the integration in l|10fl is performed over a finite interval, C = [—a, a], or 
over the complete real line, C = (—00, 00) ~ R. From (JSJ, J5J one immediately finds 

[m,cj>] v = [^,Hcj>} v . (11) 

The Krein space inner (scalar) product [t(i, <p\p has the following properties: 
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It coincides with the more general VT inner product of C. Bender el al [HI [7J 



(iPAhr = / [VT^{x)\(j){x)dx (12) 
Jccc 

when the integration path C C C of the latter integral is restricted to (an interval of) the real line, C C R, 
(i>,4>) V T= f [PTiP(x)]4>(x)dx = [ [V^(x)*}(f>(x)dx = f tt>(x)*V<Kx) = [il>,<l>]T>. (13) 

JCCS. JCCM JCCM 

In contrast to the " usual" positive definite metric structure of the Hilbert space H = TL+ © 7i _ , 

(n, (,.)), gm=/ r{x)4>{x)dx= [ (r+ct>++r-<t>-)dx, (u) 

v ' JCCR JCCR 

with non-negative norm 

" = (y^) = IKMI 2 + lhM 2 >o, (is) 



the scalar product [ip, (p]-p defines an indefinite metric structure in the Krein space IC-p = Tl + © 7i_, what 
is easily seen from the decomposition © 

[iMt>= / (r+<p+-r-<p-)dx. (i6) 

icci 

• In rough analogy with time-like, space-like, and light-like (isotropic) vectors in Minkowski space, one 
distinguishes Krein space vectors of positive type, [iJ)+,ip+]-p — ||V>+I| 2 > 0, of negative type, [ip-, ip-]j> = 
— HV'-H 2 < 0, and neutral (isotropic) vectors: 

h^b = o, v = ^+ + ^-, ||V'+II 2 = IIV'-II 2 - (17) 

In order to make the structural Krein space analogies of PTSQM models and MHD dynamo setups maximally 
transparent, we rewrite the eigenvalue problem, Hip = Eip, for the VT— symmetric Hamiltonian JSJ in an 
equivalent 2x2 matrix operator representation. Introducing the projection operators 

P±:=i(I±7>) (18) 

we decompose wave function tp and Hamiltonian H (see, e.g., [3]) as 

tp = P+tp + P-tp = 1/J+ + ip- (19) 

H = P+HP+ + P-HP+ + P+HP- + P_HP- . (20) 

In terms of the notation 

H±± := P±HP± = -d 2 x + V+(x), H ±T := P±HP T = iV-(x) (21) 

this gives 

&:)(£)-*(£)■ p =('° -')■ (22) 

where 

H±± = Hl ± , H+_ = -HL + . (23) 

If one replaces the matrix entries in (|22|l . (|23|l by appropriate constants one arrives at the schematic two-level 
model 

*»(:)-*(:)■ -V c-a)< < 2 "> 

which may be read as an elementary exemplification of Heisenberg's linear-algebraic approach |2U) to (VT— sym- 
metric) Quantum Mechanics and which was intensively studied in Refs. [51 171 IT51 |2"T1 12"2"1 l2"!i] . 
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2.2 The spherically symmetric MHD a 2 — dynamo 

The magnetic fields of planets, stars and galaxies are maintained by homogeneous dynamo effects, which can be 
successfully described within magnetohydrodynamics (MHD). One of the simplest dynamos is the spherically 
symmetric mean-field a 2 — dynamo in its kinematic regime. This dynamo model is capable to play a similar 
paradigmatic role in MHD dynamo theory like the harmonic oscillator in Quantum Mechanics (QM). Its operator 
matrix has the form 1 

d 'M = ( gfi} 1 -S[i]) ^ 



and consists of formally selfadjoint blocks 



Q{ a ]:= pap + a l l±ll, (26) 

where p = — i(d r + 1/r) denotes the radial momentum operator. The operator if; [a] is defined over an interval 
fi = [0, 1] 3 r and acts on two-component vectors <j> which describe the coupled I— modes of the poloidal and 
toroidal magnetic field components of a mean-field dynamo model with helical turbulence function (a— profile) 
air). 

Although the dynamo model is not VT— symmetric, its operator if; [a] shares a basic underlying symmetry 
with PTSQM Hamiltonians — a Z2— graded pseudo-Hermiticity ( J— pseudo-Hermiticity) [SJ^J which is induced 
by the fundamental (canonical) symmetry: 

Hi[a] = JH}[a]J, J= ( ° J V (27) 

Similar to the reflection operator V in Q , the operator J is unitary and involutive 

J t = J" 1 , J 2 = I. (28) 

The boundary conditions on the vector function <f> are set at r = 1 (the rescaled surface radius of the star or 
planet whose fluid/plasma motion maintains the dynamo effect) and it is assumed that a(r > 1) = 0. In the 
case of physically idealized boundary conditions at r = 1 (see, e.g., Ref. HH), the domain V(Hi[a\) of the 
operator -ff;[a] consists of functions <j> such that 

1 I • j, r a/ - 1/ m u u — r - /T> ^2, 



T>(Hi[a]) := U> = I j£ J : <j> eH = H®H, H = L 2 (Cl, r 2 dr), 

n=[0,l], 0(1) -0, r<Hr)U ^0}, (29) 

and Hi [a] is self-adjoint in a Krein space 

(}Cj,[.,.]j), \^,4>]J= I ^Jcj>r 2 dr, (30) 



(i 



[Hi X ,<P]j = [x,Hl<P]j. (31) 
It should be noted that for physically realistic boundary conditions 



B,. 



r=l 



= 0, B l= diag[9 r + + l)/r, 1] (32) 



there exists no appropriate Krein space which could make the operator Hi [a] J— self-adjoint. 

The structures of PTSQM models and the a 2 — dynamo can be compared most explicitly after passing from 
K.j to an equivalent Krein space with diagonal metric operator fi and redefined Hilbert spaces components, 
L2(£l,r 2 dr) 1— > £2(^5 dr). The diagonalization yields 

J~M:=($ _° I )=S-iJS, S=±(\ "J), (33) 



1 See Appendix^Jfor a few comments on the origin of this operator matrix and on the physics of a 2 — dynamos. 
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H tH » H t [a] - S ->H,W = \ ( a ) > (") 



_ 1 f 02 + ^1 



(35) 



whereas the unitary mapping [/ : L2(fl,r 2 dr) i— > L2(f2, dr) simplifies the structure of Q[a] and leads in l|34|). 
(I35[l to the additional replacements 

01.2 i-> A, 2 := r^i, 2 , Q[a] i-» <?[a] := rQ[a]r _1 = -d r a(r)d r + a(r) ^-^ ^ . (36) 

By inspection of l|21|) - (|23[1 and IM-il) - (|36|l we find that, in the chosen Krein space representations of the PTSQM 
model and the a 2 — dynamo, the block structures of the metrics (involution operators) V and fi coincide, V = ji, 
but that the blocks of the VT— symmetric Hamiltonian and the dynamo operator show significant structural 
differences 2 : 

H±± = -dl + V + {x) — _ 9[ i] ± «!«Lt£ , 

H ±T =W^ X ) ~ (37) 

It is clear that these differences in the differential expressions (as well as the different boundary conditions on 
the two-component eigenfunctions) will lead to different global behaviors of the corresponding operator spectra. 
Nevertheless, both types of systems share the same Krein-space induced features of level crossings, what will 
be briefly sketched in the next subsection. 



2.3 Spectral phase transitions 

Since the first PTSQM paper of Bender and Boettchcr it is known that VT— symmetric Hamiltonians have 
a real spectrum when VT— symmetry is an exact symmetry and not spontaneously broken 3 (the corresponding 
eigenfunctions are invariant under a VT— transformation) , whereas spontaneously broken VT— symmetry is 
connected with complex energies. A consistent PTSQM applicability and interpretation of the complex-energy 
states remains an open question up to now (cf., e.g., |26|). For convenience, we shall call these states "unphysical" 
here. 

The mathematically most interesting questions of PTSQM concern the transition between the physical 
and unphysical domains of their parameters. In the simple two-state model (I24|l it is easy to deduce that 
the quantized energies E are real ("physical") for \a\ > \b\ while they form complex-conjugate pairs in the 
"unphysical" regime where \a\ < \b\ The boundary of its PTSQM applicability coincides with the double- 

cone hypersurface in parameter space where \a\ — \b\. One easily verifies that whenever |6| approaches \a\, the 
separate eigen-energies E± as well as the corresponding two independent bound-state eigenvectors coalesce and 
coincide. On the critical hypersurface the remaining (geometrical) eigen-vector becomes supplemented by a so 
called associated vector (algebraic eigen-vector) ^3] and the Hamiltonian matrix H acquires a Jordan-block 
canonical structure ^] |57| ■ The latter cannot be diagonalized and it only gives the doubly degenerate and real 
single "exceptional-point" eigenvalue E = Ei EP ^ = c (cf., e.g., Ref. j^Hl for more details). 

An exhaustive and consistent bound-state interpretation of the Schrodinger type equation Q is more dif- 
ficult. For example, it requires the restriction of the range of exponents to a finite interval of v G (—1,2) for 
ip(x) £ L 2 {— oo,oo) as usual defined on the real line p. A rigorous proof of the reality of the energies turned 
out unexpectedly difficult ^3 EH]. For larger exponents v, the real line must be replaced by an appropriately 
deformed contour in the complex plane ^ |2] . 

A systematic analytical study of phase transition points is still lacking for PTSQM models; the same concerns 
efficient mathematical tools for deriving their location in parameter space. Similar to the double-cone hyper- 
surfaces for the simple matrix model (|24[) . one expects more complicated (and more interesting) global phase- 
transition hyper-surfaces in case of the Schrodinger type systems. Knowing the location of these phase transition 
hyper-surfaces, one would know the boundaries of the "physical" regions of exact VT— symmetry. 

2 In a very rough a nalo gy, the alpha profile a(r) has some similarities to a position depending mass, as it was studied for QM 
models, e.g., in Refs. 1251 . 

3 We recall that this follows from JU, the eigenvalue equation Hip = Eip and its VT— transformed, HVTip = E'VTtp. For 
real eigenvalues, E = E* , it is natural to set ip = VTip, whereas E ^ E* necessarily implies ip j£ VTij). 
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Figure 1: Real and imaginary components of the a 2 — dynamo spectrum as functions of the scale factor C of an 
a— profile a(r) = C x (1 — 26.09 x r 2 + 53.64 x r 3 — 28.22 x r 4 ) in the case of angular mode number I = 1 and 
physically realistic boundary conditions ll-'il'l) . The concrete coefficients in the quartic polynomial a(r) have their 
origin in numerical simulations of the field reversal dynamics (see Ref. |3(Jp. Only the imaginary components 
with 3 A > are shown. The complex conjugate (3 A < 0)— components are omitted for sake of brevity. 



For the a 2 —dynamo both types of eigenvalues — real ones as well as pair- wise complex conjugate ones 
— have a clear physical meaning. They simply correspond to non-oscillatory and oscillatory dynamo states, 
respectively. But again it is of utmost interest to know the parameter configurations for which transitions 
between the two types of states (phases) occur. In the recent paper [3U], strong numerical indications were 
presented that magnetic field reversals (interchanges of North and South poles as they are evident from paleo- 
magnetic data on the Earth magnetic field |31p are induced by a special type of nonlinear dynamics in the 
vicinity of spectral phase transition points. 

The qualitative features of the real-to-complex phase transitions are essentially the same for PTSQM models 
and for the MHD a 2 — dynamo. They correspond to transitions from Krein space states of positive and negative 
type to pair-wise neutral (isotropic) states |1 II I14| — and a square-root branching of the spectral Riemann 
surface |331 134) . Such transitions are a generic feature of Krein-space setups and they are new compared with 
setups in Hilbert spaces with purely positive metric structures as in "usual" QM. The square-root branching 
behavior is easily seen by passing from the linear eigenvalue problems for the 2 x 2— operator matrices H and 
Hi[a] of Eqs. (22) and 

(H — E)ip = 0, (Hi [a] - A) 4> = 0, (38) 

4 In the concrete case, the nonlinear transition mechanism between kinematic and saturated dynamo regime (a brief outli ne o f 
the corresponding physics can be found in Appendix [XJ was simulated with the help of a so called a— quenching (see, e.g., \32\ ) 
which simulates the nonlinear back-reaction of the induced magnetic fields on the a— profile a(r). 
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via substitutions 



^={-T^[H + + + -E}^ + )< * " ( I W (i) + A] 



(39) 



to the equivalent quadratic operator pencils 



(ff— - E)j±-(H ++ - E) - = 0, 

(Q[l] + A)-(Q[1] + A)-Q[a])0i = 0. (40) 
a J 

Both pencils are of the same generic operator type 

L[\]ip = [A 2 X 2 + AiA + Aq] tp = (41) 

with a scalar product 

(tp,L[X]ip) = a 2 A 2 +aiA + a = 0, Oj := (V>, Ajip) (42) 
which can be used to deduce the local square-root branching behavior of the spectrum 



Ai, 2 = ^ y-ai ± \Ja\- Aa a 2 j . (43) 

A typical a 2 — dynamo spectrum with a large number of real-to-complex transitions is presented in Fig. ^ 
(see also Refs. ^J[35]). These crossings with real-to-complex transition occur at exceptional points (in the 
sense of Kato 36 ) of (square root) branching type (371 155] and the corresponding eigenvalues have geometric 
multiplicity one and algebraic multiplicity two |14| . In contrast, crossings without real-to-complex transitions 
are of the same type as level crossings in Hermitian systems |11| — with geometric and algebraic multiplicity 
two |39j. Finally, we note that although locally crossings with real-to-complex transitions occur, in general, 
only between two spectral branches, globally much more branches are involved in mutual crossings (see Fig. 
^) . This reflects the fact that in general the spectrum forms a multi-sheet Riemann surface over the parameter 
space of the theory (see e.g. [3*51 l3~4l FTO] ) . 

In the next section, we will analyze the spectral behavior of a VT— symmetric interpolation model where we 
will find a similar rich structure of real-to-complcx transitions as for the a 2 — dynamo. 



3 VT— symmetric interpolation between square well and harmonic 
oscillator 

In Schrodinger-type models over the open real line x £ (—oo, oo) a PTSQM-related separation of the 
"physical" and "unphysical" domains is, in general, a mathematically highly non-trivial problem. Its resolution 
requires a fairly subtle and rigorous mathematical argumentation |lf>U29| . A typical result of the WKB analysis 
of Ref. pQ was that in a half-open interval of v £ [0, 2) the energies remain real and that the VT— symmetry 
of the wave functions remains unbroken. In parallel, a characteristic unphysical behavior of the system has 
been found in the half-open interval of v S [—1,0), where at any v < all the sufficiently high- lying energies 
E n with n > no(v) "decay" in complex-conjugate pairs, QE n ^ 0. Moreover, the spectrum becomes empty in 
the Herbst-Hamiltonian limit of the leftmost v = — 1 £Q. 

3.1 Toy model VT— symmetric differential equation 

In the present section, we are going to extend the consideration of the Schrodinger-type system JQ) to exponents 
from the interval v G [—2,0]. The end-points of this interval correspond to the purely real- valued Hermitian- 
system spectra of a freely moving particle with shifted off-set energy (for v = —2) and a harmonic oscillator 
(for v = 0). For the exponents v S (—2, 0) we expect a phase of spontaneously broken VT— symmetry with an 
involved picture of real-to-complex spectral phase transitions. 

In order to keep the numerical analysis sufficiently simple and robust, we assume the system located in a 
square well 5 (box) of finite width 2b < oo and Dirichlct boundary conditions imposed at the walls, ip(x = ±6) = 

5 Various aspects of square-well-related PTSQM setups have been earlier considered, e.g., in Refs. I11II41||42I . 
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0, i.e. we introduce an IR cut-off at the low-energy end of the spectrum. This enables us to rescale Eq. |QJ to 
the equivalent equation 

[-d 2 y + Gy 2 {iyT] i/>[x{y)] = (i(E) ^[x{y)] (44) 
with parameter-independent boundary conditions 

#K±1)]=0, (45) 

but rescaled coupling constant and energy 

G = gb A+v , ^(E) = b 2 E. (46) 

In this notation, the original bound-state problem Q with asymptotic Dirichlet boundary conditions at x — > ±oo 
is replaced by the equivalent new problem defined within a fixed finite interval [—1,1]. In the limit of very small 
b rs the potential term becomes negligible, G ~ 0, and the interaction degenerates to an infinitely deep 
square well (box) at all v. A completely similar situation occurs for systems with any non-vanishing finite 
6, but very small exponents, v rs — 2. In both extremal cases the problem remains exactly solvable. The 
original Bender-Boettcher problem corresponds to the strong-coupling limit, G — > oo, b — > oo, with g hold fixed, 
g = Gb- 4 - v = const. 

For finite coupling constants < G < oo one expects the energy spectrum to be divided into three sectors: 
into a low-energy sector with states which are involved in real-to-complex phase transitions, into an inter- 
mediate sector, where the ^—dependent energies still remain real, and into a high-energy sector with almost 
v — independent purely real eigenvalues whose states experience only a small perturbations from the complex 
interaction term. The division into low energy and intermediate-and-high energy sectors has been qualitatively 
described in a recent paper jll| by Langer and Tretter who considered a square well model with an arbitrary 
VT— symmetric potential V as perturbation. Starting from the energy spectrum of the empty square well, 
Hk = fc 2 7r 2 /4, k = 1,2,... , they showed that there are no real-to-complex phase transitions for levels k > k s 
with k s as the lowest level satisfying the supremum bound ||V||oo < (2fc s + l)7r 2 /8. In case of our model with 
potential V(y) — gb A+ y 2 (iy) v the supremum norm (see, e.g., @S]) reads (for v > —2) 

II^IU = sup \V{y)\ = \V{±1)\ = |# 4+ " (47) 
we[-i,i] 



so that it is ensured that there are no phase transitions for levels 

1 



k > k s (b) > 



4+1/ 



(48) 



According to JT] it holds for the corresponding real eigenvalues fik- (/•*& — fc 2 7r 2 /4| < |<?|6 4+I/ . The supremum 
bound is safe, but at the same time rather rough |44| . The subsequent exact numerical analysis shows that, 
depending on the concrete exponents v, the real-to-complex phase transitions in the model l|44|) stop at much 
lower energy levels. 



3.2 The emergence of SsE / on certain finite subintervals of v £ (—2, 0) 

In the generic case with v e (—2,0) and b > 0, we have solved Eqs. I)44|) + 145|) numerically by means of 
a shooting technique with a fifth-order Runge-Kutta method, utilizing and adapting standard routines from 
Numerical Recipes |45| . The corresponding code had been validated extensively in earlier work by comparison 
with known analytical results and other numerical results in dynamo theory and quantum mechanics. 

A sample of the results of such a study is depicted in Figures[21-El where we have chosen g = 1 and displayed 
the first few energy levels E{v) over the entire interval v € (—2, 0) for the sequence of values b — 2, 4, 6, 7. The 
important results of this numerical experiment are the following: 

• At all sufficiently small b, as sampled in Fig. [2 the energy spectrum exhibits a more or less ^—independent 
square- well form. 

• At not too large b the spectrum, as sampled in Fig. [31 proves clearly separated into the high-lying part 
(where the energies still preserve their approximate ^—independence), an intermediate perturbative part 
(where the perceivably ^—dependent energies still remain all real) and the low- lying part (where one 
encounters the first real-to-complex phase transitions). 
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Figure 2: Spectrum of the "PT— symmetric interpolation Hamiltonian H = —d 2 + x 2 (ix) as function of the 
exponent v for the cut-off length 6=2. All eigenvalues are real and almost independent of v. The spectrum is 
only slightly deviating from that of an empty square well. 



• The actual lowest level numbers k c (b) (critical level numbers) of the modes which are not involved in 
real-to-complex transitions lie much below the safe supremum bounds k s (b) of inequality l|48l) . Choosing, 
for example, the exponents v = —1/2 and v = —3/2 we read off that 
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• Starting from the "intermediate width" region, sampled at b — 6 in Fig. we find that the left-hand 
half of the picture exhibits a clear transition from the slightly non-Hermitian square well regime (with the 
higher energies all real) to its more strongly non-Hermitian extension where for all exponents v not too 
distant from v = — 2 the purely imaginary and finite component of the potential resembles the spatially 
antisymmetric part of the exactly solvable VT— symmetric Heavyside step potential within a square well 
considered in Ref. This explains why in Figs. |31-[S]the continuing decrease of v makes the respective 
two or three lowest pairs of the energies merge and complexify. 

• At "sufficiently large" cut-offs b, all the real low-lying energies depicted in the right-hand halves of Figs. 
0] and obviously stabilize and approach the b — > oo limiting pattern as published in Ref. [Tl 116). In 
particular, we see that the ground-state energy remains real and that it starts growing more quickly only 
when the values of v move down and closer to the Herbst limit of v — > — 1 + . We observe that in a more 
appropriate way this growing real branch should be interpreted as a special type of ladder-shaped merger of 
intermediate real segments which actually correspond to level-pairs with higher mode numbers. A zoomed 
view on this peculiarity is presented in Fig. where it is clearly visible that a chain of exceptional points 
is located on this branch with alternating complex-valued segments branching off to the left and to the 
right. These segments fit, after further complex-to-real transitions, to the real eigenvalues of the v — > 
and v — ► —2 limit models. 

• When the cut-off b is increased the following simultaneous changes in the spectrum can be observed. In 
the upper low-energy region with v > — 1 step by step more and more level pairs become twisted into the 
complex sector. With a "slight delay in 6" and at v < — 1 the lower of the twisted levels undergo a second 
pairwise real-to-complex transition with the levels below them. A sort of web structure is forming with 
a purely real branch remaining between the left [y < —1) and right (y > —1) purely complex (twisted) 
spectral regions. The complex- valued level pairs are branching off from the real branch forming a ladder- 
shaped structure. At the low-energy end of this ladder a second process occurs. The left (y < — 1) 
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Figure 3: Real and imaginary components of the spectrum in the case of a cut-off length b — 4 (complex conju- 
gate (OA < 0)— components omitted, as well as further higher lying levels without real-to-complex transitions). 
The low-energy sector with its multiple real-to-complex transitions starts to form. 



complex level pairs are passing the line v = — 1 and move to the right of it. There, at some v > — 1 
the corresponding exceptional point merges with an exceptional point of a right branch. As result one 
of the real segments between left and right off-branching levels disappears and a smooth complex- valued 
branch forms which extends over a large v interval and whose imaginary components are increasing very 
fast when b is increased. It remains the real branch which becomes more and more vertical whereas the 
complex branch is not intersecting with it (the real component of the complex branch is coinciding at one 
point with the real branch but the imaginary components are not coinciding). 

Analyzing the sequence of Figs. [2 -El we observe that, when the cut-off b is increased, a rather special (and 
seemingly inextricable) branch pattern 6 of real and complex eigenvalues is forming in the vicinity of the exponent 
v = — 1. The extreme steepness of an increasing number of imaginary branches and their accumulation at v = — 1 
(see Figs. 0JEJl as well as the occurrence of the almost vertical real branch (Fig. EJ) are indicating the formation 
of a local spectral singularity with d v E\ v ^_ 1 — > ±oo at the (almost) vertical segments of the real- valued branch 
as well as on the imaginary branches close to the exceptional points of the "ladder" structure. From the figures 
it is not at all obvious how this pattern is compatible with the Herbst limit of an empty spectrum for b — ► oo 
at v = — 1. We will resolve this interesting puzzle in the next section. 

6 The phenomenon may be generic since in Ref. 1461 . the "wiggly upwards" spectral pattern has also been detected for a very 
different one-parametric family of asymptotically exponential VT— symmetric potentials V(x) = — (isinhx)' 3 near the Herbst-like 
exponent (3 = 1. 
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Figure 4: Spectrum for a cut-off length 6 = 6. 

4 The Herbst limit and its relation to the Squire equation of hydro- 
dynamics 

In Ref. |15| it was was shown by I. Herbst that the spectrum of a Hamiltonian (|50|) with imaginary linear 
potential (imaginary homogeneous electric field) over the real line x G R is empty. The differential expression 
of the corresponding operator coincides with that of the VT— symmetric Schrodinger type equation l|44|) with 
exponent v = — 1, 

[-d*-igb 3 y]i>(y) = b 2 Eil;(y), = ±1) = 0. (50) 

The only difference of l|50ll compared to the Herbst model is in the Dirichlet boundary conditions at y = ±1 
which restrict the system to a box (square well) . Due to this analogy and for sake of brevity, we will call the 
model JH0I) a "Herbst box". 

We start our consideration by noticing that Eq. H5()(l and the spectral function /i(6, E) = b 2 E are invariant 
under the rescaling b i— ► g~ x t z b, E i— > g 2 ^ 3 E so that henceforth we can set g — 1, without loss of generality. 
The corresponding Herbst box Hamiltonian we denote as 

H Hb (b)~-d 2 -tb 3 y. (51) 

Equation l|50(l itself is of Airy type and its solutions can be expressed as 

</%) = CiAxKfoJJ + CaAaKfo)], (52) 
^(y) := e^i-iby-E], (53) 

where Ci,2 = const, and Ai(£), A 2 (£) are any two of the Airy functions Ai(£), Ai(g£), Ai(g 2 £) with q := e i27r / 3 . 
As usual, the boundary conditions lead to a characteristic determinant which defines the spectrum of the 
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Figure 5: At a cut-off length 6 = 7 the generic structure of the spectrum is clearly visible. The web-like pattern 
of the real components (Fig. Ep,) contains purely real branches in the vicinity of the left (y = —2) and right 
(y = 0) end points of the considered interval as well as a chain of purely real intermediate segments in the 
vicinity of v = — 1. (See Fig. |S|for a detailed view.) The increasing number of imaginary components (Fig. \^jp) 
with high gradients \d v E(y « — 1)| 3> 1 which accumulate in the vicinity of v = — 1 are first indications of the 
formation of a local spectral singularity at v = — 1. 



eigenvalue problem. In case of Eq. (|5U1) , it reads 

A(E)=A 1 (t + )A 2 (^)-A 1 (Z_)A 2 (£ + ) = 0, e ± :=e( y =±l). (54) 

Characteristic determinants of this type (built over Airy functions) have been intensively studied since 1995 in 
a paper series of Stepin |471 14*5] and Shkalikov et al 49 7 on the spectral properties of the Squire equation of 
hydrodynamics 8 

H Sq {e) :=isd 2 y +y, (H Sq -\) X = 0, x{v = ±1) = 0, e := {aR)~ l , (55) 

Before we make the (obviously existing) relation of this model to the Herbst box model explicit, we briefly 
review a few of its properties. 

The Squire equation l|55(l describes the normal vorticity of a plane Couette flow with linear velocity profile. 
The parameter a > denotes a real-valued wave number which originates from the decomposition of a 2D-flow- 
perturbation 

*(x,y,t) = x(y)e i&{x - xt \ (56) 

R > is the Reynolds number and e — the viscosity. The spectrum of Hs q was found to have a Y— shaped form 
All the eigenvalues are located in a close vicinity of the three segments (1, — i/V3\, ( — 1, — i/\/3], 

7 For related work see also Ref. 1501 . 

8 The corresponding physical background can be found, e.g., in 15111.521 . 
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Figure 6: The purely real curve (high-lighted fat) in the vicinity of v = —1 (cut-off length b — 7) is formed 
by intermediate real segments between complex valued segments which branch off to the left and to the right. 
This leads to a ladder-like structure with exceptional points as nodes. The zoomed picture in Fig. 03 shows 
that the "oscillations" of the real curve about the line v = — 1 persist also at its lower end, but with strongly 
reduced " amplitude" . When the cut-off is slightly increased to b > 7 the lowest exceptional (real-to-complex 
transition) point (Fig. \fijp) will cross the line v = — 1 and will in the region v > — 1 coalesce with the other 
(nearest) exceptional point. As result one intermediate segment will be removed from the real curve and a 
purely complex-valued branch-pair will smoothly tend from "far left" (v < —1) to "far right" (v > —1) 
similar to the lower lying purely complex branches visible in the graphics. 




[— i/ V3, —ioo). In the limit of large R — > oo, a > 1 and correspondingly small e — » + the eigenvalue problem 
(|55|l turns into a singular perturbation problem and its eigenvalues show a remarkable limiting behavior: For 

e — > + , more and more eigenvalues "move in" from —ioo along the line [— i/v3> — ioo), merge pairwise in 
the vicinity of the point — i/v3 and depart then (again pairwise) to move symmetrically along the segments 
[1, — i/v3], [— 1, — iA/3] and to "fill" them step by step — leaving the Y shape-invariant. The process was 
described in Ref. |48j as a special type of transition from a discrete spectrum to a continuous one. Explicitly, 
the following asymptotic estimates were found in 1^] 

A n ~ —is — - — G [— i/Vi, — ioo), tl — > oo (57) 

A± - ±l±e 1 / 3 s n e ±l7T / (i € (±l,-i/V3], e -> 0+ (58) 

where s„ are the zeros of the Airy function 

Ai(s„) = 0, s„eM_. (59) 

One clearly sees that the smaller e is chosen the smaller the distances between the eigenvalues become — leading 
in the limit e — * to a quasi-continuous spectrum. 

Noticing that the pairwise merging and splitting (level crossing) of the eigenvalues occurs at A + ~ A - ~ 
— i/v3, it is easy to estimate that the value e n , for which this crossing is connected with the n— th Airy function 
root s n , is given by 

el/ 3 \-= . (60) 

|s„|V3 

Let us utilize the above results now for the Herbst box model. A simple comparison of the eigenvalue 
problems (|50|l (for g = 1) and l|55|l shows that these problems may be made coinciding if one sets 

is' 1 [H Sq (e) - A] X (y) = V [H Hb (b) - b 2 E] V [Vr/>(y)] - (61) 
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and identifies 

b 3 =e-\ E = ib\, VH Hh {b)V = ie- 1 H Sq {e), V^{y) = X {v) ■ (62) 

This means that the two models are related by the combined action of a rescaling, a Wick-rotation and a 
coordinate reflection V . 

With the help of the estimates (|57|l . I|58|) it is now an easy task to explain the behavior of the Herbst-box 
spectrum E{b). 

• The rescaled spectrum E(b)/b = i\(e — b~ 3 ) (shown in Fig. UJ) is simply the Wick-rotated version of the 
original shape-invariant "Y" of the Squire operator Hs q (e). With increasing b — e" 1 / 3 more and more 



0.6 
0.4 
0.2 

!±l 
E 

-0.2 



X 




b=4 
b=5 


X 
X 


m 




b=6 


□ 


□ 




b=7 


o 










1 


J OEM O □ OX H 




OXCX- 










□ 








x 

















0.5 1 1.5 2 



Re(E/b) 



Figure 7: The rescaled Herbst-box spectrum E/b = i\(e — b 3 ) coincides with the Wick-rotated Y— shaped 
spectrum of the Squire operator Hs q (e). 

eigenvalues "move in" from +oo and "fill" the two complex conjugate branches (i, l/y/S\, (—i, 1/V3] of 
the "Y" as well as the half-line [l/-\/3, +oo) — in a similar way as in the original A(e — > 0) limit. For 
b — * oo the spectrum becomes quasi-continuous on the rotated "Y". 

• Due to the shape invariance of E(b)/b, the spectrum E(b) — ib\(b~ 3 ) itself inflates when b increases. It is 

located in the close vicinity of the segments (ib, 6/V3], (— ib, 6/V3], [b/y/3, +oo) and moves with b — > oo 
to infinity — leaving (for sufficiently high b) an empty region at any fixed finite distance from the origin 
of the spectral E— plane. Hence, we find (as required) that for b — > oo the Herbst box spectrum turns into 
the empty spectrum of the original Herbst model over the real line R. 

• For finite b, the asymptotic estimates (|57|) . (|58|) map into 

7T 2 fc 2 

E k - -tpt e 6x [l/V3,+oo), fc^oo, (63) 
E± - ±ib ± is n e ±i7r/6 £ bx(±i,l/V3], b -► oo , (64) 
and we identify 163|l as the pure square well spectrum 

/*(&, B*) = b 2 E k ~ n 2 k 2 /A E b 3 x [1/^3, +oo) (65) 

of the high-energy sector which is almost not affected by the PT— symmetric interaction. In contrast, 
the low-energy sector described by ()64|l shows a purely linear scaling behavior of the imaginary energy 
components Sr-E^, whereas the real components remain asymptotically fixed when b increases. This 
situation is also clearly visible from the numerical results presented in Fig. [S] The graphics of the 
equivalent spectrum ^i(b,E) 7 depicted in Fig. |5| provides a complementary description and shows how 
(for increasing b) the b— independent eigenvalues fi(b, E k ) of Eq. (|65ll leave the high-energy sector, obtain 
an explicit b— dependence in the intermediate-energy sector and finally coalesce and split into complex 
conjugate pairs. 
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Figure 8: Real and imaginary components of the Herbst-box spectrum as functions of the cut-off length b 
(complex conjugate (5A < 0)— components omitted). The asymptotical behavior of the complex- valued branches 
is clearly visible (constant real components and linear b— dependence of the imaginary components). 



• From the form of the spectral branches on the E— plane (rotated "Y") it is clear that the level-crossings 
in the vicinity of E ~ b/^/3 correspond to the typical real-to-complex phase transitions of VT— symmetric 
models in Krein spaces. With the help of relation the cut-off-scales b n and positions E n of the level 
crossings can be roughly estimated as 

fori ^ I Sn |V3/2, E n ~\a n \/2. (66) 

One can use the explicit values of these b n , (&i w 2.02, 62 ~ 3.54, 63 s» 4.78, 64 « 5.88, 65 w 6.87, be i=a 
7.81,...), to roughly derive the number of the lowest uncrossed modes in the cases b = 2,4,6,7. For 
b = 2,4,6 the result exactly coincides with the level crossing pattern shown (at v = — 1) in Figs. |2]-0J 
whereas the value b§ is clearly smaller than the actual transition value 6 5 ( c ) for which according to Figs. 
EJd and El holds b 5[c) > 7. 

For completeness, we note that the asymptotic approximation 48 of the Airy function roots 

2/3 

+ 0(n~ 4 / 3 ), 71 -> 00 (67) 

together with (|66|) yields the following rough estimate for the lowest purely real-valued mode 

4/2 \ 3/2 1 

fca> ^b 6 ) + 2- (68) 
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Figure 9: The rescaled Herbst-box spectrum fj,(b, E) = b 2 E allows a complementary view on the transition from 
the high-energy sector to the intermediate and low-energy sector. 

The scaling dimension n a = 3/2 of this bound is only one half of the scaling dimension n s = 4 + v = 3 of 
the corresponding supremum bound (|48(l . 

The exact positions of the level crossing points are given by the multiple roots of the characteristic deter- 
minant, A(.E) = 0, 8eA(E) = 0. This equation system can be simplified via Wronskian VF[Ai(.), ^4 2 (.)] 
to yield the conditions 

A 1)2 (£+)=±j1 1)2 (£_) (69) 

(see 0Hl f° r the details). 

Plugging the numerical results from the eigenvalue solver into this equation with j4i j2 (£) chosen as in 
Ai(£) := Ai(£), A 2 (0 := Ai(g 2 £), selects the condition Ai i2 (£+) = Ai, 2 (£_) and satisfies it within 
numerical working precision. For the same data holds Ai j2 (£ + ) 7^ — Ai. 2 (£_). 

• The spectral behavior E(b) for increasing cut-off b can be summarized as follows. At the beginning, the 
real eigenvalues from the high-energy sector decrease as Ek ~ it 2 k 2 / (4b 2 ) — moving into the intermediate 
energy region. When b approaches b n ~ |s„|\/3/2 from below, the real eigenvalues {i? 2n _i, -E 2n } (cor- 
responding to a pair of positive and negative Krein space states |11| 1 coalesce at Em-\ ~ -E 2n ~ l s n|/2 
and a real-to-complex transition occurs {£J 2n _i, S 2n } — ► , When 6 is further increased the real 

energy components remain fixed WIE^ w |s n |/2 (see Fig. EJi), whereas the imaginary components blow up 
linearly along the asymptotes ~ ±&± s n V3/2 (Fig. Eb). 

Let us, for finite b, relate the obtained Herbst-box results to the spectral behavior of the VT— symmetric 
interpolation model of the previous section. Apart from the obvious one-to-one correspondence of the high- 
energy sectors (see (JHEJ), a clear identification is immediately possible for those Herbst-box eigenvalues which are 
close to the imaginary axis and which have the largest imaginary components. These eigenvalues are located on 
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the branches with the largest imaginary components in Figs. (which stay complex when v passes through 

the Hcrbst-box value v = —1). It is clearly visible from these figures that, for increasing b, the imaginary 
components are blowing up, whereas the real components remain asymptotically constant. 




V 



Figure 10: Concrete example for the generic merging process of two complex-valued spectral branches (present 
here for b < b coa i « 6.36) into a single complex- valued branch for b > b coa i. The two exceptional (real-to-complex 
transition) points existing for b < b coa i are located on the plane (v^b^E — 0) and coalesce at b — b coa i. The 
corresponding point (v coa i ~ —0.9983 > — 1, b coa i 6.36, ^sE = 0) is the starting (cusp) point of a sharp "valley" 
of the non- vanishing imaginary component which steeply grows and smooths when b > b coa i is further increased. 
The white curve marks the Herbst-box values. Before the two exceptional points coalesce, the left one of these 
points crosses the Herbst-box configuration at (y = —1, b = b c sa 6.02, *<sE = 0) and a real-to-complex transition 
occurs for the Herbst-box model. 

So far, we have found a clear correspondence for those regions on the Herbst-box "Y" which are located 
away from the center of the "Y" with its real-to-complex phase transitions. A more subtle situation occurs 
in the vicinity of this center. The corresponding Herbst box eigenvalues will map into points located close to 
(or on) the forming (almost) vertical segment of the purely real branch depicted in Fig. From the zoomed 
graphics in Fig. we observe that the purely real and almost vertical branch of the interpolation model 
"oscillates" around the Herbst box line at v = — 1 with strongly decreasing "amplitude" to its low-energy 
part. With increasing b this "decreasing amplitude" effect becomes stronger and the "oscillations" are only 
traceable with the help of an appropriately increased zooming scale. Nevertheless, it can be read off that the 
real-to-complex transition of the Herbst box spectrum follows qualitatively the same scheme for any finite b. 
The real eigenvalues of the Herbst box are all located on the purely real branch of the interpolation model 
and the real-to-complex transition occurs when the lowest exceptional point on this branch moves from the left 
sector v < — 1 through the Herbst-box value v = — 1 into the right sector v > — 1 — to coalesce afterwards 
with the next higher exceptional point from the right sector. With this passing of the left-sector exceptional 
point through the line v — — 1 the Herbst box eigenvalues become pairwise complex conjugate with strongly 
increasing imaginary components (due to the asymptotically diverging gradient d v E{b — > go)\ v ~_y — > ±oo). 
The real-to-complex transition with subsequently increasing imaginary components are illustrated in Fig. 1101 

Finally, we note that in the limit b — > oo, the lowest-lying intersection of the purely real branch with the 
Hcrbst-box line v = — 1 moves away to infinity like b/^/3 (the lower bound of the real segment [b/V3,oo) of 
the Herbst-box "Y") so that the real branch itself remains for any finite energy in the right sector v > —1 — 
approaching the Herbst-box line asymptotically. This reproduces the earlier observations of Refs. P Ej for 
the spectrum of the Bender-Boettcher problem over the real line. Additionally, our Herbst box results predict 
for this problem diverging imaginary components at v = —1: |3_E(6 — > co,v — > — 1)| —* oo. Taking these 
observations together we once more see that in the limit b — ► oo a spectral singularity is forming at v = — 1 
with \E(b — > oo, v — > —1)| — > oo, \d v E(b — > oo, v — > — 1)| — > oo. 
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5 Conlusions 



In the present paper we considered three models emerging in different physical setups, but which are closely 
related with each other by their underlying mathematical structure as spectral problems in Krcin spaces. The 
models are a one-dimensional VT— symmetric quantum mechanical interpolation setup defined over a square well 
of finite width 2b, the spherically symmetric MHD a 2 — dynamo as well as the Squire equation of hydrodynamics. 
For the PTSQM model and the a 2 — dynamo we made their close relation transparent by transforming them 
into a 2 x 2 matrix operator representation with coinciding block structure of the Krein space metric (involution 
operator). In the case of the Squire equation we showed that the corresponding spectral problem is connected 
with a VT — symmetric eigenvalue problem by a rescaling and Wick rotation 9 . 

Based on recent results on the spectrum of the Squire equation, we performed a qualitative analysis of the 
VT— symmetric quantum mechanical interpolation model for arbitrary square well widths (cut-offs) 2b. This 
allowed us to trace the emergence of the Herbst limit with its empty spectrum as a spectral singularity and 
to fit our results to those of the Bender-Boettcher equation over the real line. We obtained a rich structure of 
multiple spectral phase transitions from purely real eigenvalues to pairs of complex conjugate ones — as it was 
to expect for spectral problems in Krein spaces. 

A deeper insight into the Herbst-box spectrum and a possible extension of the present results to VT— sym- 
metric Hamiltonians of the type Hm,n = —d% + x 2M (ix) N , M, N = 1, 2, 3, . . . over square wells can probably 
be achieved by representing the characteristic determinant A(E) in (|54|l via Hadamar product representation 
of the Airy functions |53| as a spectral determinant of Bethe-ansatz type 1161 1291 1541 and studying it by similar 
cocycle functional equations as in Ref. |55| . 

A question which was not touched in the present paper concerns the orthogonality of the Herbst-box eigen- 
functions. For the Squire equation it is known that its cigcnfunctions show a strong non-orthogonality | 47l 15 1| 
(due to the non-normality of the Squire operator) for eigenvalues in the vicinity of the branch point center of 
the "Y" (pseudo-spectral techniques HOI IS EH play an important role in this case). Our above considerations 
indicate on a link of this issue with the forming spectral singularity \d v E(y —1,6 3> 1)| S> 1 in the vicinity 
of the almost vertical segments of the purely real branch in the spectrum of the VT— symmetric interpolation 
model. 

Finally, we would like to note two issues which seem of relevance for future considerations. The first one is 
in developing efficient mathematical tools to find the hyper-surfaces in parameter space where spectral phase 
transitions of the real-to-complex type occur 10 . Knowing these hyper-surfaces, one would know the boundaries 
which separate the parameter space regions with unbroken VT— symmetry from regions with spontaneously 
broken VT— symmetry. In case of a 2 — dynamos the corresponding knowledge would allow for a more precise 
prediction of configurations with tendency to magnetic field reversals. The second issue concerns methods 
for solving inverse spectral problems in Krein spaces. Such methods would be extremely helpful for the data 
analysis of the dynamo experiments which are planned for the near future at seven sites around the world |58| . 
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A A few comments on the physics of MHD a 2 — dynamos 

The dynamo operator Hi[a] originates from the MHD mean-field induction equation (cf. \12\ ) 

d t B = V X («B) + v m AB (70) 

for the magnetic field B. This equation results from averaging over small scale turbulences in the velocity field 
of the electrically conducting fluid (or plasma) which drives the dynamo. The helical turbulence function a(x) 
(also called a— profile 11 ) encodes the net effect of the small scale physics on the large scale (mean) magnetic 
field B. For certain topologically non-trivial helical velocity and B— field configurations an inverse cascade 

9 It is clear that, apart from the Squire equation, there will exist other hydrodynamic equations which can be structurally 
identified as Wick-rotated VT— symmetric systems in Krein spaces. 

10 A two-step method similar in spirit was successfully used, e.g., in higher-dimensional gravitational models to obtain the 
stability regions in the moduli (parameter) space of these models (step on e: fi nd the critical hyper-surfaces; step two: identify the 
stability/instability properties of the model aside of these hyper-surfaces) 1571 . 

11 In general setups, a is not a scalar function but a tensor 1121 . 
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effect occurs which induces an energy transfer from small-scale structures to large-scale structures (inverse 
to the energy transfer in usual turbulence cascades where the energy is pumped from large-scale structures 
into smaller structures until it finally dissipates and transforms into thermal energy). For sufficiently strong 
inverse cascade effects the advection term Vx(aB) starts to dominate over the diffusion term ^ m AB (y m is 
the magnetic diffusivity) and the magnetic field strength starts to grow exponentially. This kinematic dynamo 
effect (growing B— field for a given velocity field of the fluid) is followed by a saturated dynamo regime where 
a balance between the dynamo effect and the back-reaction of the induced magnetic field on the velocity field 
(via Navier-Stokes equation) prevents a further growth of the field strength B. For completeness, we note that 
an MHD dynamo is an open system in which part of the kinetic energy of the conducting fluid (or plasma) 
transforms into magnetic field energy. 
The dynamo eigenvalue problem 

Hi[a]<f>i >n = Xi, n (pi, n , fa,n(t) ~ expA ; ,„ t (71) 

follows from the induction equation 1)70(1 via a double decomposition: decomposing the B— field into poloidal and 
toroidal components (what leads to the two-component vector structure of 4>i,n) and expanding them further 
into spherical harmonics. In the simplest (toy model) case of a spherically symmetric dynamo configuration the 
corresponding modes decouple completely and one arrives at the spherical I— mode projection l|25|l. (|71() (the 
subscript n denotes the radial mode number). 

Up to now only a single exactly solvable a 2 — dynamo model is known — the model with constant a— profile 
[T2| . Its spectrum is discrete, real [211, bounded above and, depending on the value of a, it is either completely 
negative (for a below a critical a c : a < a c ) or it contains a finite number of positive eigenvalues A; in > 0. 
The dynamo effect is dominated by these latter eigenmodes. In practice, it usually suffices to concentrate the 
analysis on the dominating upper most growing mode (or a few of the upper most modes) for dipole (1=1) and 
quadrupolc (1=2) configurations. (There exist no "s— wave" a 2 — dynamos with I = 1121 I13| .) 

The spectral properties of the dynamo operator Hi [a] are becoming much richer for inhomogeneous a— profiles 
a(r) ^ const, when real-to-complex transitions occur — as discussed in Subsection 12.31 and shown in Fig. ^ 
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